Implementation of quantum lògic operations 
and creat ion of entanglement 
in a silicon-based quantum computer with constant interaction 

G. P. Berman 1 , G. W. Brown 2 , M. E. Hawley 2 , D. I. Kamenev 1 , V. I. Tsifrinovich 3 

1 Theoretícal Division T-13, Los Alamos National Laboratory, 
■ Los Alamos, New Mexico 8754-5, USA 

o 

çS| ■ 2 Materials Science and Technology Division MST-8, 

\ Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA and 



(N 



(N 



3 



IDS Department, Polytechnic University, Brooklyn, New York 11201, USA 



We describe how to implement quantum lògic operations in a silicon-based quan- 



tum computer with phosphorus àtoms serving as qubits. The information is stored 

"xl" ■ 

in the states of nuclear spins and the conditional lògic operations are implemented 



through the electron spins using nuclear-electron hyperfine and electron-electron ex- 



. change interactions. The electrons in our computer should stay coherent only during 



implementation of one Control-Not gate. The exchange interaction is constant and 
selective excitations are provided by a magnètic field gr adient. The quantum lògic 
operations are implemented by rectangular radio-frequency pulses. This architecture 
is scalable and does not require manufacturing nanoscale electrònic gates. As shown 
in this paper parameters of a quantum protocol can be derived analytically even 



cr. 

> 

for a computer with a large number of qubits using our perturbation approach. We 
present the protocol for initialization of the nuclear spins and the protocol for cre- 
ation of entanglement. All analytical results are tested numerically using a two-qubit 
system. 

PACS numbers: 03.67.Lx, 75.10.Jm 



I. INTRODUCTION 



The long decoherence time of nuclear spins of phosphorus donors in silicon makes quantum 
computers based on these spins attractive for quantum information processing. A scanning 
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FIG. 1: A schematic illustration of two phoshorus àtoms placed in a permanent magnètic field B°. 
The electron spins S\ and S2 (large arrows) of the neigboring àtoms interact with each other via 
the exchange interaction and the nuclear spins I\ and I2 (small arrows) interact with the electron 
spins through the hyperfine interactions. 

tunneling microscopy technique can be used for creation of many identical 

arrays of phosphorus àtoms on the (100) surface of silicon. The phosphorus qubits can be 
encapsulated 
time. Kane 



3V overgrowing additional silicon layers 5] to increase the electron relaxation 
6[ proposed to use nanoscale electrònic gates to control the qubits. This 
technique has not yet been realized, due to fabrication issues, and so in this paper we 
consider a different architecture. In our approach, the exchange interaction between qubits 
is constant and selective interactions are realized through the use of a magnètic field gradient 
and both microwave and radio frequency pulses. Measurement can be implemented using 
optical techniques . The measurement can be facilitated by creation of many identical, 

noninteracting spin chains to amplify the signal. 

The 31 P atom has electron spin 1/2 and nuclear spin 1/2. If the qubits in each chain are 
placed at a separation of ~ 10 nm from each other, the nuclear-nuclear, nuclear-electron, 
and electron-electron dipole-dipole interactions are small compared to the electron-electron 
exchange interaction so that one can neglect the dipole-dipole interactions (see Fig. EJ). 
There is also a relatively strong hyperfine interaction between the electron and nuclear spins 
of a 31 P atom. Since the relaxation time for the electron spins at temperatures of 1-7 K 
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is relatively small (0.6-60 ms at 7 K the quantum information must be stored in the 

states of the nuclear spins. Because the nuclear spins do not interact, electron spins can be 
used to mediate the nuclear-nuclear interactions. In this setup, the electron spins must be 
coherent only during the relatively short time of implementation of a quantum lògic gate, 
such as a Control-Not gate, on a particular pair of qubits. 

We consider in this paper a procedure for implementation of entanglement between the 
nuclear spins in a two-qubit quantum computer using rectangular radio-frequency pulses. 
Entanglement is the simplest quantum lògic operation needed for implementation of more 
complex quantum lògic gates and is useful for demonstration of quantum computation in a 
potentially scalable solid-state system. The paper is organized as follows. The eigenstates 
of the system are calculated analytically in Section II. A brief description of the protocol 
for creation of entanglement is given in Section III. The quantum dynamics of the system 
is described in Section IV. The eigenstates from Section II are used for calculation of pulse 
parameters in Section V. Initialization and entanglement with two qubits are simulated 
numerically in Section VI. In Section VII we review the working conditions and the parameter 
range for our computer. 



II. EIGENSTATES 



The unperturbed Hamiltonian of the system reads 



H° = le B{S{ + le B z S z - ln B{l{ - ln B z 2 I z 2 + 



A 



1 



sm + s z 2 r 2 + - (s+iï + sïi+ + s+i 2 + s 2 i+ 



j 



QZ qz , 1 ( q+ q— , Q— Çf-f 
°1°2 2 \ 1 2 1 2 



where S% (1%) is the projection of the kth electron (nuclear) spin on the ath axis, k — 1, 2, 
a = x,y,z; S^ 1 = ± iS%, if = ± iï\] 7 e and 7 n are, respectively, the electron and 
nuclear gyromagnetic ratios; B^ is the permanent magnètic field in the location of the fcth 
spin; A and J are, respectively, the hyperfme and exchange interaction constants. 
Let us define 

± = si + s z + i z + i z . 

The quantum states form 5 independent subspaces characterized by the value of 



S n = (if) n \T,\i/j n ). 
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Two one-dimensional subspaces with S n = ±2 are formed by the eigenstates with the 
following eigenvalues E n and eigenvectors \i/j n ): 



#15 = \ {-leBt - J e B* + ln Bl + ln B Z 2 ) + - + j = (- 7e + ln ) b + - + 



E = (7e 



|« = | | 2 «1 II), 
ln)b+^ + j, |-0O> = 



T2t2^ltl), 



(1) 

(2) 
(3) 



where b = (_B| + -Bf)/2 and we use the notation | fï^) and | Jjfc) for the states of the electron 
spins and | an d I for the states of the nuclear spins, k = 1, 2. 
The basis vectors for the subspace with S n = — 1 are 



(4) 









1 12^2^1 |l) 




|H) 


| Í2«l tl) 




|14) 


J T2JJ-2 JJ-1 ll ) j 




{\D j 



Introducing the notation SB 
( 



[B 2 ~ the Hamiltonian matrix for E„ 



-j e óB + 7„ò - 



A 
2 







J 

2 
2 



2 



V 



le 5B + ln b-{ 

-7e& + InàB + I 



•1 becomes 
\ 

■ (5) 



f - 7e 6 - 7n<^ + f y 

This matrix can be diagonalized analytically. For 5B = the eigenvalues and eigenfunctions 



were found in Ref. ll|. Instead of the exact analytical solution we apply here a perturbative 
approach. Our perturbative approach has an advantage over the exact analytical solution 
because it can be used to find the eigenstates for the quantum computer with more than two 
qubits when no exact analytical solution is available. The perturbation theory is based on the 
fact that 7 e 6 is at least three orders of magnitude larger than J/4, A/2, ^JSB, ^ n b, and ^ n 5B. 
(We take b w 3.3 T so that 7 e 6/(27r) » 93 GHz, J/(2vr) «1-10 MHz or less, A/(2n) = 116 
MHz.) For our range of parameters the matrix (J3J) splits into two relatively independent 
2x2 blocks . The first block is formed by the matrix elements in the upper left còrner 

and the second block is formed in the lower right còrner. The relative independence of the 
two different blocks follows from the facts that (a) the moduli of the differences between 
the eigenvalues of each block are much smaller than the moduli of the differences (~ 7 e 6) 
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between the eigenvalues of the different blocks; and (b) the matrix elements A/2 relating 
the different blocks are much smaller than 7 e 6. The corrections to the wave function due to 
the neglected terms are of the order of 

(A/2) _ 1A _ 4 



e = 



6 x 10" 



(6) 



and corrections to the energy levels Ef', i = 1, . . . , 14, are of the order of tens of 
kilohertz. These corrections are important in order to flip the nuclear spins because the 
Rabi freqüències of the nuclear spins are of the same magnitude. 
The eigenvalues and the eigenfunctions are 



,(o) 



E. 



(0) 



E { $ = - le b + ln 5B + \^>) = | | 2 h ^ U >, 



J 



l^ 0) ) 



■,(0)\ 



E 



(°) - n, h - - 

13 — 7n° ^ 



E 



(0) 

íi 



J 

4 

J 



j 2 



lnb-J + d(l e 5By+ . 



The eigenfunctions corresponding to E$ and e[°J are 



le 5B + J( le 5B)* + Ç 



(7) 
(8) 

(9) 

(10) 

(11) 



J, 



le 5B + J( le 5By + Ç 



fr^i)>®Ui), (12) 



where i — 1, . . . , 14, are the normalization constants. The correction E^ is calculated 
in Appendix A and all corrections Ef\ i — 1, 14, are listed in Appendix B. 
The basis vectors for the subspace with E n = 1 are 







(\2)\ 


1 T2^2^ltl) 




|4) 


1 Í2M1M 




|1) 


J |2^2fTl Tl > y 




{\S}) 



The Hamiltonian matrix for E n = 1 has the following form: 



<y e SB -inb-{ 



j 

2 

A 
2 



--f e 5B - 7n& - T 



A 
2 





A 
2 



A 
2 



7e ò - ln 8B + { 




le b + ln 5B + {) 



(13) 
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The eigenvalues and the eigenfunctions are 



E<®= le b- ln 5B + j, \rj>\ 



J 



(0)\ 



t 2 t 2 ^i li >, 



Ef ] = 7e ò + ln 5B + \^>) = I 12 ft 2 ÍTltl), 



J 



,(0)\ 



|^ 0) ) = I>2| Í2 



7e 55+ J( le 5B)i + Ç 



l4 0) ) 



7^ + ^/(7^ + ^ 



(14) 
(15) 
(16) 
(17) 

lTi>, (18) 
2 1Ti> > (8) I Ti>- (19) 



The six-dimensional multiplet with £„ = splits into three relatively independent sub- 
spaces. The first two eigenvalues and eigenfunctions are 



^6 



- le b - ln b + |4 0) ) = I Í2 ^2 h Tl ), 



(o) 



(20) 
(21) 



The remaining four eigenstates with the basis vectors 







(m) 


1 Í2^ 2 ^lll) 




|5> 


1 12^2^1 Tl) 




|10) 


v | Ï2M1 Tl) j 




U 2 )J 



are related to the following Hamiltonian matrix: 

j 
2 

<5 J B(-7e-7n)-4-4 



/«5B( 7e -7n) + 4-Í í 



V 



J 
2 











2 4 

SB^e +7 n; - 2" - 4 

4 







J 

2 



,J5(-7 e + 7n) + f -Í / 

(22) 



One can see that these states form two independent (in our approximation) two-dimensional 
subspaces. The first subspace is described by the 2x2 matrix in the upper left còrner and the 
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second subspace is described by the 2x2 matrix in the lower right còrner of the matrix 
The eigenvalues and eigenfunctions are 



E 



(0) 



J í A\ P 

- Jn SB-- + \ h e SB + -) +—. 



(23) 



|4 0) )=^3|Í2)® 



.4 



A\ 2 P 



le 5B + - + \ [ le 5B + -) + — 



J 



\^ ] )=D,\h) 



E 



J 



(o) 



J 



4 



A\ 2 P 



ln 8B---\[ le 8B + -\ + — 



.4 



A\ 2 P 



le 5B + - + \ / 7 e íS + - + — 



(0) 



J 



.10 - Tri&B - - - y (7e&B - - 



|V>S?> = £>ioiI 2 > 



le 5B-^ + \l[ le 8B- 



J 



^ ( 3) = D 12 \ | 2 ) 



1)2 h ) + -I U-2 iT .1. ) } : ! Ii), 

(24) 
(25) 

|^1h)J®Ui>, 
(26) 
(27) 

|«i)U|íi>, 

(28) 
(29) 

J 



A\ 2 P 



A\ 2 P 
+ — 



A\ 2 P 



1\2 h ) + -| hl\l) }®\W- 

(30) 

The eigenvalues E\o and E\2 in Eqs. (|27j) and ()29jl are written for the case A/2 > 7 e <5-B. For 
the opposite case, A/2 < jeóB, one must exchange the eigenvalues E w <-> i? 12 and leave the 
eigenvectors l^io) an d \1p12) unchanged. 



III. CREATION OF ENTANGLEMENT 



Consider the quantum dynamics generated by electromagnètic pulses for different param- 
eters J and óB. The four states, 6th, 7th, 14th, and 15th, have the lowest energies of the 
order of — 7 e 6. The distance between the lower 6th and the upper 15th levels of the quartet 
is (2 7n 6 + A)/(2?r) « 173 MHz [7„/(27r) « 17.25144 MHz/T]. This is much smaller than 
/cbT = 20.83 GHz, where T = 1 K is the temperature. Consequently, all four of these states 
are initially populated. 
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The initialization of the nuclear spins and creation of entanglement between the nuclear 
spins can be implemented by using the fact that the electron spins are polarized. Our system 
can be represented as a one-dimensional spin chain 

|n 2 e 2 eini). (31) 

In Eq. (JÏÏTJ) ei assumes the vàlues ÏTi , Jji , e 2 =fÏ2 , -JJ-2 , ™i =Ti , li , and n 2 =] 2 , 12 , In the 
spin chain (|31|) there are interactions only between the neighboring spins, so that this kind 
of spin ordering is convenient for analysis of conditional quantum lògic gates. 

Initially our chain is in the superposition of states with the lowest energies Eq, E7, En, 
and E15. From Eqs. (j2J), (|ZJ), ©, and (ffilj) one can see that these states are 

\n 2 hh ni), (32) 

with different n\ and n 2 . One possible initial state is shown in Fig. ^ By using Control- 
Not gates between the electron and nuclear spins, one can transfer the polarization from 
electron to nuclear spins. After some time the electron spins polarize and one obtains the 
only populated state 

I | 2 «i!i>. 

By using the Hadamard transform on the lst nuclear spin, Control-Not gate between the 
lst nuclear spin and lst electron spin, Control-Not gate between the lst electron spin and 
2nd electron spin, and Control-Not gate between the 2nd electron spin and 2nd nuclear spin 
one can create entanglement between all spins of the system 

-^(| UÍl2ÍliU)+e ie \ Uhhii))- (33) 
The exact value of the phase 9 is not important for us. 

IV. DYNAMICS 

The time- dependent magnètic field has the following components: 

B 1 ^) = B l (cos(vt + <p),- sin(z/í + <p),0), (34) 

where B l , v, and ip are, respectively, amplitude, frequency and phase of the pulse and t is 
time. The frequency v can assume both positive and negative vàlues as shown in Fig. [21 
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FIG. 2: Different kinds of polarization of the electromagnètic wave fbr v > and v < 0. 



The perturbation term in the Hamiltonian has the form 



V(t) 



S 1 + S 2 



e -i(yt+<p) _|_ 



(35) 



where = 7 e E 1 and Q° = 7„5 1 . 



Scheme for numerical simulations 



The numerical simulations are performed without using the perturbation approach and 
results are presented in Sec. IVIl below. It is convenient to work in the rotating frame where 
the effective Hamiltonian is independent of time. The relationship between the wave function 
in the laboratory frame and the wave function in the rotating frame is 



!*(*)> 



!*(*)>• 



The Schròdinger equation in the rotating frame is 



(36) 



z|$(t)) = (H° + vZ + V) |$(í)), (37) 

v = n e (sï + si) - nS(i? + íj). (38) 

Let us decompose the wave function over the eigenstates \ip n ) of the Hamiltonian Hq as 

|$(í)> = £c„(í)|^n>, (39) 



n=0 
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where the functions \ip n ) are related to the basis functions |i) by 

IVO = Ç an.il*)- (40) 

i 

The coeficients a n ^ are calculated in Section [iï] in the zeroth order approximation. (For 
our numerical simulations we use the exact vàlues of a n ^.) The system of 16 differential 
equations for the expansion coeficients c n (t) is 

15 

ic n (t) = {E n + E n u) c n (t) + V nm c m (t), (41) 

m=0 

where 

V n , m =(i/> n \V\i/> m ). (42) 
Eq. (|41jl can be regarded as the Schròdinger equation 

15 

^n(^) ^ ] h"n,mCm 

(t) (43) 

m=0 

with the time-independent Hamiltonian h whose matrix elements have the following form: 

The dynamics of the coeficients c n {t) can be computed using the eigenfunctions b q n , q = 
0, . . . , 15 and the eigenvalues e q of the Hamiltonian h as 

c n (t) = £ c m (t ) £ (blY ble-^l (45) 

m=0 q=0 

where t' is the time of the beginning of the pulse. The wave function in the laboratory frame 
can be represented as 

\^{t)) = Y.D P {t)e-^\%). (46) 

Before each pulse at time t = t' we make the transformation to the rotating frame using the 
relation [see Eqs. I|BE |) . ([SÜ Jl and (@ïïj)] 

D p (t) = exp [iE p t + %L v {yt + ip)\ c p (í), (47) 

and after the pulse at the time t = t' + r (r is the duration of the pulse) we make the back 
transformation to the laboratory frame using the same formula. 
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V. IMPLEMENTATION OF LOGIC GATES 

We describe here implementat ion of lògic gates in terms of the basis states \i), i = 
0, . . . , 15. From Section |H] the basis functions \í) are approximately equal to the eigenfunc- 
tions \tpi) if the conditions jeSB 3> J/2 and \^ e àB — A/2\ 3> J/2 are satisfied. Assume that 
the frequency of a pulse is close to the transition frequency of the kth (k = 1, 2) electron or 
nuclear spin; the direction of the spin in the state \p) is along the direction of the permanent 
magnètic field B° (i.e. \p) = \ . . . . . .) or |p) = | . . . f k . . .)); and the state \q) is related to 

the state \p) by a flip of the kth. spin. For the initial conditions 

D p (t') = 1, D q {t') = 

the dynamics of this spin is described by the following equations : 

-iA qiP r/2 



|cos 


Aç,pT 


. -^q,p ■ 

+ 'í^- L·í - sm 


\q,pT 


j 




2 




2 





D q (t' + r) 



n 



x 



sm 



q.p 



^q,p T 



iA qtP t' -iip+iA q , p r /2 



(48) 



where 



^q,p Eq Ep ^q,P 



i! is the time of the beginning of the pulse, r is the duration of the pulse, Q = Q e is the 
Rabi frequency of the electron spin, and Q = Q n is the Rabi frequency of the nuclear spin. 
For the other initial conditions 

D p (f) = 0, D q {t') = 1 



the solution is 



D p {t' + 



D q (t' + T ) = { cos 





= i- 


— sm 




|cos 


\,pT 


2 





\,pT 



i- sm 



A 



Ag,pT 



D iA qtP r/2 



(49) 



The complete transition between the states \p) and \q) takes place when the detuning is 
equal to A q>p = and when r = tt/Q (7r-pulse). A near-resonant transition with A^O can 
be completely suppressed when the condition jïj] 

IA„ 



y/AK 2 - 1 ' 



(50) 
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known as the 2irK condition, is satisfied. Here K = 1, 2, . . . is an integer number. For this 
value of fl the value of the sine in Eqs. (jüj) and (jlüjl is equal to zero. 

As follows from the above considerations, (i) in order to implement a complete transition, 
the frequency of the pulse must be resonant to this transition, v = E q — E p , where v can 
assume both positive and negative vàlues, (ii) In order to completely suppress a transition 
with A ^ 0, the Rabi frequency of the pulse must satisfy the 27rií-condition (J3ÜJ). Both 
operations (i) and (ii) can be implemented simultaneously by one pulse if there are two states 
in the quantum register. Actually, we described here the procedure for implement ation of 
the Control-Not gate, which we will use to create entanglement in our system. 

We now derive the parameters of the gates described in Section IIHI Let CNj^ be the 
Control-Not gate which flips the kth spin in the state \i) and suppresses the flip of the same 
spin in the state \j) (the latter has different orientation of the Zth. spin). We assume that 
the kth spin is pointed up (along the direction of B°) in both states \i) and Let the 
state \i') be related to the state |z) by the flip of the kth spin and the state \f) be related 
to the state \j) by the flip of the same spin. Then the frequency v and the detuning A^/j in 
Eq. (JïïUJ) for the Rabi frequency are 

v = Ei> — Ei, Aj/j = Eji — Ej — v. (51) 

In our notation, it is convenient to treat the state of the electron spin | ff-^) as |0fc) and the 
state | JJ-fc) as The fact that the energy of the state |0fc) is larger than the energy of the 
state |lfc) is accounted for by a negative frequency v of the pulse in the first equation (|5Tj) . 



A. Initialization 



The nuclear spins can be polarized by using the fact that the Larmor freqüències of 
the electron spins depend on orientations of the corresponding nuclear spins through the 
hyperfine interaction. Measuring the electron Larmor freqüències using, for example, a 
scanning tunneling microscope 16, Q], one can define the orientation of the nuclear spins 
and apply a selective 7r-pulse if necessary. Here we describe a different technique. We assume 
that initially the nuclear spins are not polarized, the electron spins are polarized, and there 



are four states (|3Zj) in the register. If we apply the gates 

CN ei , ni CN niiei (52) 
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(the order of implementation of the operators is from the right to the left), we obtain the 
superposition of states 

|n 2 JJ-2 ei |i ) (53) 

with indefinite orientation of the first electron spin and with definite orientation of the first 
nuclear spin. 

Since the state of the first electron spin in Eq. is unknown, we cannot immediately 
swap the states of the second electron and nuclear spins because of the interaction between 
the electron spins. One has to wait while the electron spin polarizes again (during, for 
example, the time-interval 0.1 s). In our numerical simulations the relaxation of the electron 
spins is modeled by flipping them "by hand", without using electromagnètic pulses, in all 
states of superposition. 

After the electron spins are polarized one applies the gates 

CN e2in2 CN n2 , e2 (54) 

and obtains the state 

\Ue 2 hli)- (55) 
One waits while the second electron spin relaxes and obtains the state 

|| 2 ^i|i> (56) 
which is used as an initial state for creation of the entanglement. 



B. Entanglement 

The sequence of gates 

CN e2 , n2 CN ei , e2 CN m , ei Had ni (57) 

generates the entangled state (|33|). In Eq. (|57j) Had ni is the Hadamard gate on the first 
nuclear spin, and the gate CN^y is the inverse of the Control-Not gate: it flips the target 
qubit only if the control qubit is in the state |0fc). The parameters of these gates can be 
calculated analytically using Eq. ([51jh In our simulations presented below, the Hadamard 
gate is performed by applying a n/2 pulse of duration r = 7r/(20 n ), where Q„ is given by 
Eq. (JEU) below. 
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In principle, it is possible to inplement a Control-Not gate between the nuclear spins 
without changing the states of the electron spins In practice, this approach is not 

useful because one can show (using calculated in this paper eigenvalues) that the mediated 
by the electrons effective coupling between the nuclear spins is of the order of 3.5 Hz or 
less. This means that the Rabi frequency of the pulse implementing the Control-Not gate 
must be less than 2 Hz and the frequency of the pulse must be tuned in resonance with the 
accuracy of approximately 0.1 Hz. 



C. Rabi freqüències 

The Rabi frequency of the electron spins Q e is different from Q° = ^ e B l because of 
the exchange interaction between the electron spins. Similarly, the Rabi frequency of the 
nuclear spins O n is different from 0° = jnB 1 because of the hyperfme interaction between 
the nuclear and electron spins. 

Consider, for example, the transition 

\hhhh)^\hhhh), (58) 

associated with the flip of the first electron spin. The matrix element = í) e /2 of the 
matrix h njTn in Eq. (|44|) is responsible for this transition. The value of íl e must satisfy the 
2tïK condition (|50Jl . The value of Vq^ can be calculated using Eq. (}4~2*j) . Only two terms in 
the sum 

15 

V nm = a i,n(Í\V\j)a m} j, (59) 

i,j=0 

where n = 6, m = 4, appreciably contribute. These are 

V 6A = Y o»,6<i|V'|j)o4^ ~ a6,6(6|V"|4)a 4i4 + a 6 ,6(6|V"|2)a 452 . (60) 

i,j=0 

The first term is due to the flip of the first electron spin by the electromagnètic pulse, i.e., 
due to the transition 

| Í2«l íl)^|Í2^líl>, (61) 

The second term is due to the two-step transition 

| Í2«líl> - I Uh^lU) - I T2^2^ltl), (62) 
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where the first step is due to the nonselective excitation of the second electron spin and 
the second step is implemented due to the exchange interaction between the electrons. 
The contribution of the second term is proportional to the ratio of the matrix element 
J/2, responsible for the exchange interaction, [see the second transition in Eq. Í)fi2j1 ] to the 
detuning 25B, characterizing the first transition. From Eq. (|HÜ|) the Rabi frequency is 

^K 1 -^)- < 63) 

This expression for the electron Rabi frequency holds for all other electron transitions. The 
2i\K condition (|5üj) for the electron spin becomes 

flQ = 7e fli = l A ^l (l L\ ( 64 ) 



where |Aj/j| « A when the control spin is the nuclear spin and |Aj/j| « J when the control 
spin is the neighboring electron spin. 

We now find the Rabi frequency of the nuclear spin. Consider the gate CN ei m in Eq. ([52)1 
acting on the state | ^2 JJ-2 -U-i Ti )• The nuclear spin is flipped as a result of the transition 

\hhhu)^\hh^iïi), (65) 

which is implemented through the matrix element Vqj = —íl n /2. In Eq. (|42j) . only two 
terms considerably contribute to the value of Vqj, namely 

V 6 j « a 6i6 ( — ^ I a 7J + a 5fi ^-a 7J . (66) 

The first term is due to the flip of the nuclear spin by the electromagnètic pulse in Eq. í)65j) . 
The second term is due to the two-step transition 

| Í2«líl> - I 12^2^1 |l) ^ | hhhh), (67) 

where the first step is due to the hyperfine interaction between the first electron spin and the 
first nuclear spin and the second step is due to the nonresonant action of the electromagnètic 
pulse on the first electron spin. The transition ()67|) is initiated by the electron spin, which 
creates a magnètic field in the x direction comparable to the magnètic field B 1 of the pulse. 
In spite of the fact that the probability P e of flipping the electron spin is small, the fact that 
Vl e /VL° n ~ 10 3 makes the probability of the transition (J67|) comparable to the probability of 
the transition (|S5|). 
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In order to calculate a 5 6 in Eq. (jSSjl. we must calculate the first-order correction to the 
wave function (J2~üj) using the equation 
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— P (0) p (0) 



Taking n = 5 and multiplying by (6| we obtain 



15 V n 

F (0) F (o) ■ ( 69 ) 

771^5 

Only one term with m = 6 significantly contributes to the sum, where 

^5,6 « (5|f |6) = |. 

We find 

05,6 * 2[4 0) - E®] * 2 7e 6 <<C L 
Putting this value to Eq. (jüüjl we obtain the expression for the Rabi frequency of the nucleus 

n n = n° n + ^- b n° e « 2í£. (70) 

This equation is vàlid for the Rabi freqüències associated with other nuclear spin transitions. 
Similar to Eq. (jfiHj) for the electron spin, the correction to the Rabi frequency is proportional 
to the ratio of the matrix element A/2, responsible for the hyperfme interaction, to the 
detuning 7 e 6 « 7 e 6 — 7 n 6 between the freqüències of the nuclear and electron spins. 
The 2nK condition for the nuclear spin reads 

^ = InB' = )M (l + ^rY\ (71) 



where \^f,j\ ~ A. We have calculated the Rabi freqüències only for the two gates. The Rabi 
freqüències for the other gates can be calculated using the same formulas (JHlj) and (fTTf . 

We still have indefinite parameters K e and K n in Eqs. (|51j) and (|71|) for the Rabi freqüèn- 
cies of electron and nuclear spins. Increasing K e and K n decrease the electron and nuclear 
Rabi freqüències which should satisfy the conditions 

«1, £n = -r%«i. (72) 



A le 5B A ln 5B 
These conditions provide selective excitations of the spins The error in the probability 
amplitude due to nonselective excitations is of the order of e e for the electron spins and of the 
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f the 
22, 



order of e n for the nuclear spins. The condition (J72j) can be written as B 1 / (4(5 B) -C 
distance between the qubits is 10 nm and the magnètic field gradient is 10 5 T/m [20 
then B 1 < 46B = 2 x 10~ 3 T. Hence, 

^ w < 56MHz, (73) 

27T 27T 

We now choose K = K e for the electron spin to satisfy the following two conditions: (i) 
The value of K e must be large enough to satisfy the first equation (JZ2J), which allows one to 
decrease the error due to nonselective excitations; (ii) and the time r e of implementat ion of 
the Control Not gate on the (target) electron spin must be much smaller than the electron 
relaxation time T c (we assume T c = 60 ms ^(|). This condition can be satisfied by decreasing 
K e . 

Consider the Control-Not gate on the electron spin with control nuclear spin. The error 
\[P~ e due to nonselective excitation of the electron spins is of the order of [for |A| fa A in 

Eq. K5DJ] 

. p — ça ça (7A\ 

V A le 5B 8 le 5BK e K e ' K J 

where we used the parameters A/2n = 116 MHz, the magnètic field gradient 10 5 T/m, and 
the distance 10 nm. In our simulations we take K e = 103 so that the error is \[P~ e ~ 0.01. 
The Rabi frequency and the time duration of the 7r-pulse are 

a A 



2?r 2itxklC* - 1 



563 kHz, r e w 0.89 fis. 



Consider the Control-Not gate on the electron spin with control electron spin. The error 
is of the order of (|A| fa J) ^j~P' e ~ 0.01 if K' e = 1, where we assume J / (27r) = 500 kHz. The 
Rabi frequency and the time duration of the 7r-pulse are 

Q' J 

— fa fa 289 kHz, r ' = 2 us. 

2tt 2 7 r v /4^ - 1 

Due to the second equation (|72|) . the Rabi frequency fl n must satisfy the condition 

^ « < 34.5kHz. (75) 

27T 27T 

For implementation of the Control Not gate, íl n must also satisfy the 2nK condition í[71|). 
The error in the Control-Not gate on the nuclear spin with control electron spin is \[P~^ ~ 
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0.05 if K = K n = 33, 620. The Rabi frequency and the time duration of the 7r-pulse acting 
on a nuclear spin are 

q A 

1.7 kHz, r n ^0.3ms. (76) 



2tt - 1 

We note here that the magnètic field gradient practically defines the clock speed of our 
quantum computer. During the time of implementat ion of the Control Not gate on the 
nuclear (target) spin, the electron (control) spin must stay coherent so that the condition 
T n <C T c , where T c is the transverse relaxation time, must be satisfied. For isotopically 
purified 28 Si, the relaxation time T c can be as long as 60 ms at 7 K, which is large enough 
for implementation of the Control-Not gate. In natural Si (4.7% of 29 Si) the value of T c is 
smaller than 0.6 ms at 1.6 K ^(|, which is too small for implementation of the Control-Not 
gate on the nuclear spin. The Rabi frequency of the order of 1.7 kHz is large in comparison to 
the dipole-dipole interaction between the electron and nucleus of the neighboring phosphorus 
àtoms, which is close to 32 Hz for the distance 10 nm between the qubits, so that one can 
neglect this dipole-dipole interaction. 

We now make some comments about precision of calculation of the energy levels in 
Section|H] As follows from Eqs. (|48J) and 1)490. the spin rotates around the x axis (in the 
rotating frame) with the frequency 



K P =sjK P + W- (77) 

If we wish to implement the resonant transition, A 9iP must be equal to zero. Since A q iP = 
E q — E p — v is defined by the distance between the energy levels (eigenvalues), the complete 
transition takes place only if these eigenvalues are exactly known. For a system with a small 
number of qubits, the eigenvalues can be calculated numerically with a high precision. For a 
system with a large number of qubits, (> 30) one can use the perturbation theory described 
in Section [Hi This theory allows one to calculate the eigenvalues with precision e m A/2, 
where e is defined in Eq. (JBJ) and m is the order of the perturbation theory (m = 2 in our 
paper and m = 1 for the zeroth order approximation). The limited accuracy of E q and E p 
results in a finite detuning of the order of 



\ i 

I a <?,pI ~ e y 



In order for this detuning to have a small influence on the dynamics, the value of |A giP | in 
Eq. (177)1 must be much smaller than the value of íl. From Eq. (JïïJ) we have e = 6 x 10~ 4 . In 
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the zeroth order approximation (m = 1), the value of |A çp /(27r)| is of the order of 35 kHz, 
which is much larger than Q n / (2ir) in Eq. (fTüj) and is not acceptable. In the second order 
approximation (m = 2), we have |A 9iP /(27r)| ~ 21 Hz, which is much smaller than Í2 n /(27r) 
so that this is an acceptable approximation for us. In practice, the above argument means 
that the frequency v of the electromagnètic wave must be tuned with an accuracy of the 
order of several tens of Hertz in order to flip a nuclear spin without generating a substantial 
error. 



D. Spin relaxation 

The relaxation of the electron spin can affect (flip) the nuclear spin via the hyperfine 
interaction. Here we discuss the conditions required to suppress the nuclear spin flip during 
the electron spin relaxation process. Other kinds of relaxation mechanisms can be neglected 
because the relaxation time T n of a neutron is at least four order of magnitude larger 

n 

(T n ~ 3 x 10 3 s |2í|) than the total time of implementation of the algorithm. The electron 
spin relaxation allows us to implement the nonunitary transformation 

A)| U) -> I U) 

/ (78) 
Di\U) , 

which is necessary for creation of the initial state. 

Consider the dynamics of the classical nuclear magnètic moment I(t) placed in a perma- 
nent external field B° oriented along the z axis and a field created by the electron spin S(t) . 
Here the z component S z (t) of the electron spin is a given function of time t. The slowly 
varying z component of the magnètic field acting on the nuclear spin is 

B'(t) = B°- A;S z {t). (79) 

The electron spin rotates with the frequency 7 e -B° and generates a circularly polarized 
time-dependent magnètic field acting on the nuclear spin via the hyperfine interaction. Since 
7e-B° is three orders of magnitude larger than the Larmor frequency of the nuclear spin, this 
fast field does not affect the dynamics of the nuclear spin because the nuclear spin is out of 
resonance with the field generated by the fast-rotating electron spin. 
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There are always stray magnètic fields in a real system, such as, for example, the Earth's 
field, which can affect the dynamics of the nuclear spin. Without loss of generality, we 
assume that this field, B x , is oriented along the x axis. The nuclear spin is flipped if two 
conditions are satisfied. (a) The B z field Íl79j) must go through the zero point B z (t) = 0. 
(b) The condition of adiabatic passage j^J 

B z 

Í = -^7^2 « 1 (80) 



must hold. Here £ is a dimensionless small parameter. We now analyze how to choose the 
parameters of our system in order to violate these two conditions and, thus, to suppress the 
nuclear spin flip. 

Condition (a) is satisfied if the value of B z in Eq. (J79|) is always positive, which yields 
the minimum value -B^in °f the external field (for S z = 1/2) 

B ° > B -°- = éfi * 3 - 36 T - (81) 

Even if condition (a) is not satisfied, the nuclear spin cannot flip if the magnètic field B x 
is sufficiently small to violate condition (b). With a good approximation the function S z (t) 
can be chosen in the form 2^] 

S*(t) = 1 -(l-2±). (82) 



2 V T, 



The parameter £ in Eq. (|5ü|l becomes 



,. 6.2 x 10 3 

t « - rz , 2 , (83) 
T r ( B x 



where the dimensionless relaxation time T c is equal to the number of milliseconds [T c = 
T c /(1 ms)] and the dimensionless B x field B x is equal to the number of gauss [B x = 
B x /(l gauss)]. We take T c = 6 (so that T c = 6 ms). If B x is equal to the Earth's magnètic 
field, B x = 0.5, the condition of adiabatic passage is not satisfied, £ ~ 4.1 x 10 3 ^> 1 so that 
the nuclear spin does not flip. The numerical modeling of the classical spin dynamics with 
B x = 0.5 gauss and 

B° = 3.3T < B^ a , T c = 6 ms, I z (t = 0) = -0.5 (84) 

yields 

= i»(0) - /■(!■.) 1Q _ 4 
/•0 
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One can show that the error in the quantum probability amplitude is of the order of 5I Z . 
One can neglect this error if it is small in comparison with the other errors. For B x = 5 gauss 
we have £ ~ 41, and a numerical simulation gives SI" ~ 0.075. For B x = 50 gauss we have 
£ ~ 0.41 and numerical simulations show that the nuclear spin flips and I Z (T C ) sí 0.5. This 
is the situation when the conditions of adiabatic passage are satisfied which prevent the 
initialization of our computer. 

If the transverse magnètic field is relatively strong, for example, when B x = 50 gauss, one 
can suppress the flip of the nuclear spin by violating condition (a). For example, increasing 
B° from B° = 3.3 T, which is less than B° min in Eq. (JEU), to B° = 3.5 T > B^ in and for the 
same transverse magnètic field B x = 50 gauss, we numerically obtained 5I Z (T C ) = 7 x 10~ 4 , 
i.e. the nuclear spin actually does not flip. 

Next we will show that, in the situation when the conditions of adiabatic passage are 
satisfied for certain spins, the error still can be small if B° is close to B^, i.e. when 
K-BjJdn - B °)/ B íin\ < !• Assume that we are dealing with an ensemble of identical spin 
chains as mentioned in the introduction, B° < B^ in , and condition of adiabatic passage (fKÜj) 
is satisfied. Since before the relaxation the electron spins point in a random direction, not 
all of them pass the point for which B z (t) = in Eq. í)79|l. In Fig. we show how the 
magnètic field acting on different nuclear spins changes with time. Only those nuclear spins 
flip for which the condition 



holds. As follows from the figure, the total number of such spins in the ensemble is 



The error \fP t in the probability amplitude for the ensemble is \fP r ~ i). For example, for 

B° = 3.3 T (and for -B^ in = 3.36 T) we have rj = 0.009 which is a small error. In summary, 
our analysis shows that it is possible to suppress the flip of the nuclear spin during the 
relaxation of the electron spin and to implement the transformation (|78j) . 




(85) 





VI. NUMERICAL RESULTS 



In our four-qubit system, all parameters of the pulses can be calculated numerically us- 
ing exact eigenvalues of the matrix H° and off-diagonal elements of the matrix V, which 
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B z (arbitrary units) 



B°+A/(2r n %) 




t (arbitrary units) 
► 



B»-A/(2y„h) 



FIG. 3: A schematic illustration of the magnètic field acting on the nuclear spin for different vàlues 
of S z (0) [S z (0) = 1/2 for the lowest curve and S z (0) = — 1/2 for the top curve]. The nuclear spin 
flips for the magnètic fields illustrated by the dashed curves and does not flip for the magnètic 
fields shown by the sòlid curves. 

are equal to fl e /2 for electron transitions and — Q n /2 for the nuclear transitions. In spite 
of the ability to calculate the parameters numerically, we calculate them analytically using 
our perturbative approach and simulate the quantum dynamics numerically as described 
in Sec. IIV Al Our analysis has the following advantages: (i) it can be applied to a system 
with an arbitrary number of qubits (ii) it allows one to take into consideration only "slow" 
transitions with small detunings and to neglect fast transitions with relatively large detun- 
ings, which have little influence the quantum dynamics. Using our approach it is possible 
to understand the most important sources of error and to minimize them by the optimal 
choice of pulse parameters. 
We start with the state 



|tf (0)> = C 6 (0)|6> + C 7 (0)|7> + C 14 (0)|14) + C 15 (0)|15) 



(86) 



with arbitrarily chosen complex coefficients C 6 (0), C 7 (0), C u (0), and Ci 5 (0) at time t — 0. 
Then we make the transformation to the representation of the Hamiltonian H° [see Eq. (jHS|) ] 
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fl B (o) = E«^,(o). 



(87) 



i=0 
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FIG. 4: The probability error Fasa function of J/(2ir); 5B = 0.5 mT. 

After initialization of the system and creation of entanglement, we make the back transfor- 
mation 

C l (T) = jra hn D n (T)e-* E " T , (88) 



n=0 



where T is the total time of implementation of the protocol. The error is calculated as 



P 



-\Co(T)f 



+ 



-|C 15 (T)p 



(89) 



In Fig. E]we plot the error after implementation of initialization and entanglement. Each 
point on the plot is the average over 100 realizations with randomly chosen complex co- 
eficients C 6 (0), C 7 (0), Ci 4 (0), and Ci 5 (0). One can see that the error increases with J 
increasing. When J is large enough, J ~ 47 e í5, the basis states \n) differ considerably 
from the eigenstates \ip n ) of the Hamiltonian H°, so that the error is generated as a result 
of free evolution of the basis states. The error bars are the consequence of the fact that the 
protocol processes some initial states of the superposition better than other states. 

In Fig. El we plot the probability error P as a function of the magnètic field difference 
25 B = B>2 — Bl for the interval 0.1 mT < 25 B < 5 mT. For a qubit spacing of 10 nm, this 
interval corresponds to magnètic field gradients from 10~ 4 to 5 x 10~ 5 T/m. As follows from 
the figure, the error is large when 25 B is small, i.e. when J/{A^ e 5B) ~ 1. 

The maximum in P near 25B k, 4.2 mT is defined by the condition ^ e 5B = A/2. When 
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FIG. 5: The probability error P as a function of magnètic field difference 25B for J /{2n) = 2 MHz. 
The data are averaged over 100 realizations with randomly chosen normalized complex coeficients 
C 6 (0), C 7 (0), C 14 (0), and C 1B (0). 

this condition is satisfied the lOth and 12th eigenvectors defined in Eq. (J28j) and (|3U|) become 
symmetric and antisymmetric superpositions 



| 2 > ® 



^(|^ 2 lti>±|lMi» 



® I Ti)- 



These states are formed because the electron frequency difference caused by the magnètic 
field gradient is compensated by the hyperfine interactions between the electrons and nuclei. 

The numerical results indicate that it is worthwhile to place neighboring qubits at larger 
distance d from each other. This gives one the following advantages. (i) At a given gradient, 
the value of SB = [B 2 Z — B\)/2 increases with increasing d, which provides better selectivity 
of the pulses and increases the clock speed of the quantum computer [see the second equa- 



. This does not 



tion (fT2~j) and Eq. (JTüj)]. (ii) The value of J decreases with increasing d 
affect the clock speed of the quantum computer because the clock speed is defined by 5B. 
(iii) Decreasing J decreases the influence of the off-diagonal components of the exchange 
interaction and the eigenstates \ip n ) of the Hamiltonian H° are better approximated by the 
basis states \n). Since \ip n ) ^ \n), the error in the probability amplitude [proportional to 
J I (A'~f e 5B)} is generated even in the stationary system when no electromagnètic pulses are 
applied. (iv) In a system with more than two qubits increasing d decreases unwanted effect 
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of long-range interaction between distant (not neighboring) qubits. (v) When the distance 
between the qubits is large, the error is less sensitive to the random qubit displacements 
caused by imperfect qubit positioning using a scanning tunneling microscope 27| . 



VII. SUMMARY 



We described how to implement quantum lògic operations in a silicon-based quantum 
computer with phosphorus àtoms serving as qubits. The lògic operations can be implemented 
in our computer if the following conditions are satisfied. 

1. The selective excitations of nuclear spins can be implemented if their Rabi freqüències 
are small [see the second equation (JZ2J)], 

ü n <^2 ln {B z 2 -Bl), (90) 

i.e., when B 1 <C B% — Bf, where B 1 is the amplitude of the radio- frequency field 
and £?2 — ^ s the magnètic field difference equal to the product of the magnètic field 
gradient and the distance between the qubits. Condition (JÜÜj) defines the clock speed of 
our computer. The Control Not gate between the nuclear spins is implemented during 
the time-interval approximately equal to tcn ~ Tr/^n- The time- interval required to 
flip the electron spins is at least two orders of magnitude smaller than tcn- As follows 
from Eq. (Jüüj) . the clock speed is practically defined by the magnètic field gradient: 
the larger is the gradient, the faster is the computer. 

2. The time-interval tcn must be much smaller than the electron relaxation time T2 
because the electron spins should stay coherent during implement ation of the Control 
Not gate. 

3. In order to make the electron spins polarized, the magnètic field B z must be large 
and the temperature must be small, i.e., the condition ■y e B z 3> /cb@ (&b is the 
Boltzmann constant) must be satisfied. 

4. The electron-electron exchange interaction must be small in comparison with the fre- 
quency difference 7 e (5|— B\ ) between the electrons, i.e. the condition 7 e (i?|— B\ ) J 
must be satisfied, otherwise, the off-diagonal components of the exchange interaction 
would modify considerably the basis states and generate error. 
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5. The modified inequality Y^eiB^ — Bf) — A\ ^> J, which includes the hyperfine interaction 
constant A, must hold also. 

6. In order to suppress the flip of the nuclear spins during the relaxation of the electron 
spins, the computer must be shielded from external stray transverse magnètic fields, 
B x , so that B x < 5 gauss or the external permanent magnètic field B° must be close 
to or larger than -B min = 3.36 T. 
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APPENDIX A 



Here we calculate the corrections E^ to the eigenvalues for some states using per- 
turbation theory Q|. The corrections to the other eigenvalues are calculated in a 

similar fashion. The second order correction E^ to the eigenvalue E^ is 

K^l^ (2) l^)l 2 



(2) 



E 



E, 



(0) 



(Al) 



The matrix elements of the operator are [see Eq. (JSJ)] 



\H^\J) 



( | ^ 
| 
| 



f 



(A2) 



The basis vectors \j) are defined in Eq. (j3J). 

The transformation from the eigenfunctions \ip n ) to the basis vectors \i) is described by 
the matrix a n ^ in Eq. (j40jl . We have 



(A3) 
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From equation (JJJ) we have a i 7 = 8^ where óij is the Kronecker delta-function, and Eq. (|A3|) 
becomes 

(^i^ {2) iVn) = E^(7|^ (2) U)· (A4) 
j 

From Eqs. (@J) and (|A2[) . the only nonzero matrix element is 

(7|^ 2 )|11) = ^. (A5) 



From Eqs. (jAljl . (IÀ4I) . and (IA5|) . we obtain 

(2 ) _ A |ai M i| 
^ 7 ~ 4 4°) _ #(°) ' ( j 

For our range of parameters, 1 011,11 1 2 ~ 1. Using Eqs. (jZJ), (|TÜ|) . and ()A6jl . we obtain Eq. (jB2j). 

(2) 

Next we will calculate the correction E\ 3 . 

p (2) _ V \(^3\H^n)\ 2 

13 ^ p (o) F (o) • 

The matrix elements are 

a i,13 a n,j (^ (2) |J) = 

y 

■^"(Oll, 13^71,7 + Ol3,13 a n,14) = -^-( ail > 13( W a 13,13^n,14) • 

The second-order correction is 



^13 - 4 



1 12 I 12 

I Q-11,13 | | a 13,13| 

p(o) ^(o) ^(o) 



2 



4 P (o) P (o)- 

41 — -C/1 4 



13 — ^7 -^13 ~~ -^14 J -^13 ~~ ^14 

(0) 



In order to calculate the correction to the eigenvalue E$, we note that the state |6) is 
related by the matrix elements A/2 to the states |5) and 1 10} [which can be obtained from 
the state |6) by swapping the states of ith nuclear and electron spins, i = 1, 2]. After a brief 
calculation, one can obtain Eq. ífïïTjl. 

APPENDIX B 

The corrections of the second order Ef \ i = 1, 14, to the energy levels E? ] are 



E ( ? 



A 2 1 



4 (7e + 7n)& - InàB - J{ le 5Bf + .P/A + J/2 
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E. 



(2) 



A 2 



4 "(7e + 7n)& + 7n^ + J{ le 5B) 2 + J2 /4 _ J /2 



£Í 2) « 0, 



4 2) 



A 2 



1 



4 -(7e + ln)b - iJB - ^d e 5B) 2 - J 2 /4 + J/2' 
A 2 



1 1 

+ 



üg» - Ef ] Ef ] - E^ 



[the vàlues Ef\ i = 1, . . . , 14 are defined in SectionITT]. 



EP 



A 2 

T 



i i 

+ 



4 2) 



A 2 



1 



(2) 



4 -(7e + ln)b - ln 5B - ^d e 6B) 2 + J 2 /4 + J/2 ' 
A 2 1 



4 (7e + ln)b + iJB + J{ le 5B) 2 + J74 + J/2 



i 2) 



i! 

4 

A 2 



1 1 

+ 



F (0) p (0) 1 p (0) p (0) 
Eg - E 5 E g - L· lQ 

1 1 

+ 






P (2) 
^11 



E 



(2) 
13 



(2) 
14 



A 2 


1 




4 


(7e + 7n)& + 7ní5 + yj ' { le ÒB) 2 + J 2 /4 - 


J/2' 




« 0, 




A 2 


1 




4 


(Te + Tn)& - 7n55 - ^/( 7e 55) 2 + J 2 /4 - 


J/2' 


A 2 


1 




4 . 


-(le + ln)b + ln 5B + J^eóB) 2 + J 2 /4 - 


h J/2 



(Bl) 
(B2) 



In the text we assume = £1- °' ) + E^ . 
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